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We assume that the early evolution of matter produced in relativistic 
heavy-ion collisions is described by the transverse hydrodynamics. In this 
approach only transverse degrees of freedom are thermalized, while the lon- 
gitudinal motion is described by free streaming. When the energy density of 
the system drops to a certain value, the system becomes isotropic (locally, 
in the momentum space) and freezes out. The sudden isotropization transi- 
tion is described with the help of the Landau matching conditions, and the 
freeze-out is modeled with THERMINATOR. Within this scenario one is 
able to reproduce in the quite satisfactory way the transverse-momentum 
spectra, the elliptic flow coefficient V2, and the HBT radii of pions and 
kaons studied at RHIC (Au+Au collisions at -^/snn = 200 GeV). The ob- 
tained results indicate that the system produced at RHIC does not have 
to be fully thermalized in the early stage. On the other hand, the final 
three-dimensional thermalization is necessary to reproduce the HBT ra- 
dius Along. 

PACS: 25.75.Dw, 25.75.-q, 21.65.+f, 14.40.-n 

1. Nowadays, the early evolution of matter produced in relativistic 
heavy-ion collisions is described most successfully by the perfect-fluid hydro- 
dynamics [HEJE]. Using the appropriate equation of state and incorporat- 
ing resonances in the modeling of the freeze-out process, the hydrodynamic 
approach describes very well the hadron transverse-momentum spectra and 
the elliptic flow coefficient V2- Interestingly, with a suitable modification of 
the initial conditions, the hydrodynamic approach describes also the HBT 
radii [3]. 
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In spite of its clear successes, the use of the hydrodynamics is still chal- 
lenged by the problem of early thermalization - the correct description of 
the data requires rather early starting time of hydrodynamics [5], which im- 
plicitly assumes a very early equilibration time. A possible solution to this 
problem was proposed in Ref. [6\. With the assumptions that only trans- 
verse degrees of freedom are thermalized and the longitudinal dynamics is 
represented by free streaming (a physical picture explained naturally in the 
string models one can obtain the parton transverse- momentum spec- 

tra and v 2 which are consistent with the data [§] . We shall refer to this 
approach below as to the transverse hydrodynamics^. 

In this paper we develop the ideas and formalism introduced in [6j 
9j. The transverse-hydrodynamics model is supplemented with the sud- 
den isotropization transition followed by freeze-out. This allows for a direct 
comparison of the model predictions with the hadronic data. In particular, 
besides the transverse-momentum spectra and v 2 studied already in [6j[9], 
we analyze now the femtoscopic observables. The isotropization transition 
is described by the Landau matching conditions, and the freeze-out is mod- 
eled with the help of the Monte Carlo code THERMINATOR [12J. We 
note that THERMINATOR simulates also the hadronization process, since 
it "replaces" fluid elements with the hadronic constituents. 

Certainly, our treatment of the isotropization transition and freeze-out 
is very much simplified. More elaborate approaches would describe this 
kind of transformation using kinetic theory or dissipative hydrodynamics, 
for example, see [131 ES ESj- We expect, however, that our scheme may 
be regarded as an approximation to more advanced approaches where the 
considered processes have gradual character. 

2. The equations of the transverse hydrodynamics follow from the 
energy-momentum conservation law, dJF§ = 0, with the energy-momentum 
tensor defined by the formula [16] 

if = — [(e 2 + P 2 ) U»U V - P 2 (<T + VV) ] . (1) 
r 

Here r = \/t 2 — z 2 is the longitudinal proper time and no describes the den- 
sity of transverse clusters in rapidity. The clusters are formed by groups of 
partons having the same rapidity. They are two-dimensional (2D) objects, 
whose thermodynamic properties are described by the 2D thermodynamic 
variables: e 2 , P 2 , s 2 and T 2 (2D energy density, pressure, entropy density, 
and temperature, respectively). These quantities satisfy the standard ther- 
modynamic identities: e 2 + P 2 = T 2 s 2 , de 2 = T 2 ds 2 , and dP 2 = s 2 dT 2 (here 

1 We note that our formulation differs from the formalism of the transverse hydrody- 
namics introduced originally in [101 111] . For a more detailed discussion of this issue 
see [pj- 
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the baryon chemical potential is set equal to zero). The definition of the 
energy-momentum tensor ([1]) includes also the two four-vectors, 

= («o cosh?7, u x , u y , uq sinhr/), = (sinhr/, 0, 0, cosh ry), (2) 

where = uj_,0) is the hydrodynamic flow in the central plane where 
z = 0, while rj = 1/2 In ((t + z)/(t — z)) is the spacetime rapidity. The 
four-vectors and satisfy the following normalization conditions 

U»U» = 1, W% = -1, 17>% = 0. (3) 

The four- vector combines the motion of the fluid element in a cluster 
with the motion of the cluster, thus it corresponds to the flow four- velocity 
in the standard hydrodynamics. The term V^V V in ([1]) is responsible for 
vanishing of the longitudinal pressure, i.e., in the local rest-frame of the 
fluid element, where we have = (1,0, 0, 0) and = (0, 0, 0, 1), one finds 
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(4) 



Very interestingly, such a structure of the energy-momentum tensor appears 
in the theory of the color glass condensate and glasma for r S> l/Qs, where 
Q s ~ 1 GeV is the saturation scale expected at RHIC [TTJ, [T5] . This gives 
further support for consideration of the transverse hydrodynamics as the 
appropriate description of the early evolution of matter. 

We solve the equations of the transverse hydrodynamics numerically 
using the equation of state P 2 = £2/2 = v g T% /(2vr), where v g = 16 [16] . 
Our investigations are restricted to the midrapidity region (z ~ n ~ 0), 
where the partonic system may be treated as boost-invariant. The initial 
conditions assume that the 2D initial energy density at r = t\ = 1 fm is 
proportional to the mixture of the wounded- nucleon density, pw (x±), and 
the binary-collision density, ps (x±), 

i ^ \ v g T 2 3 (rj, x±) 1 — k . ._, . 

£2(Ti,X_Lj = ^ OC piy (X±) + Atpe (xj.) ■ (5) 

TT Z 

The distributions pw (x±) and are calculated for a given central- 

ity class from the Glauber model. Following the PHOBOS studies of the 
centrality dependence of the hadron production [19] we take k = 0.14. The 
normalization constant required in ^ determines the 2D initial central 
temperature of the system, T 2 \ = T2(ri,0). 
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3. The sudden isotropization transition is described by the Landau 
matching condition 

ZfU = ?fU, (6) 

where T^ v is the standard energy-momentum tensor of the perfect-fluid 
hydrodynamics 

If = (e 3 + P 3 )U^-P 3 g^. (7) 

Here £3 and -P3 are the three-dimensional (3D) energy density and pressure 
of the system just after the isotropization transition. Equations §§§ and (J7J) 
give 

— £2 = £3, (8) 
r 

which should be supplemented by the requirement of the entropy growth in 
the isotropization transition, 

^2 < -3, (9) 

r 

where S3 is the 3D entropy density. Dividing both sides of Eqs. ((HJ) and © 
one obtains 

T2 > p-. (10) 

In our calculation, £3, P3, and S3 are interpreted as the quantities char- 
acterizing the hadron resonance gas in equilibrium at the 3D temperature 
T 3 {. This means that we identify the isotropization transition with the 
hadronization process; if the energy density of the system drops to £3^3 f), 
the transverse gluonic clusters change into a locally isotropic resonance gas. 

The transverse hydrodynamic evolution of matter produced in the non- 
central (c = 20-30%) Au+Au collisions at ^/snn = 200 GeV is shown in 
Fig. [TJ The solid lines indicate the 2D isotherms defined by the condition 
T2 = const. The numbers at the isotherms show the values of the ratio £ = 
T2/T21. We stress that for massless particles our hydrodynamic equations 
exhibit scale invariance — the temperature field may be rescaled by an 
arbitrary factor without the change of the flow. As described below, this 
property helps us to fulfill simultaneously the Landau matching conditions 
for the energy and entropy. The dashed line in Fig. [1] indicates a possible 
transition curve defined by the condition T3 = const. As expected, the line 
of constant T3 does not agree with any of the 2D isotherms. We note that 
all the curves in Fig. [1] are plotted for = (the in-plane direction) . Due 
to the initial asymmetry of the system in the transverse plane, for <f> ^ 
the isotherms have slightly different shapes. 
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Fig. 1. 2D isotherms describing purely transverse hydrodynamic expansion in 
the Minkowski space (solid lines). The initial conditions correspond to the non- 
central (c = 20-30%) Au+Au collisions at ^/snn = 200 GeV. The numbers at the 
isotherms give the values of the 2D temperature in units of the initial central tem- 
perature T%i. A possible isotropization transition curve defined by the condition 
Tsf = constant is represented by the dashed line. 

Our model has essentially three free parameters: T21, T 3 f, and no (the 
initial time t\ is always set equal to 1 fm). In the fitting procedure we first 
assume uq = 1 and try to adjust T21 and T 3 f in such a way that a good 
description of the data is achieved. In this step, the energy-conservation 
condition ([5]) is used as the only constraint. In the next step we rescale 
the 2D temperature in such a way that the equality T2 = 3^3/(253) holds 
at the initial time r = t\. In this way we make certain that the condition 
(|10p is fulfilled on the whole transition hypersurface. The rescaling of the 
temperature changes the energy density on the 2D side, however, this may 
be compensated by the appropriate change of no- In this way we are able 
to satisfy finally the conditions (|SJ) and fl§|). 

4. Our description of the isotropization process is quite extreme. Yet, 
we add to it another extreme assumption: the evolution of the isotropic 
(locally equilibrated) system is so short that the isotropization process may 
be followed immediately by freeze-out. Technically, this assumption is re- 
alized by the Monte-Carlo simulations done with THERMINATOR, where 
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Fig. 2. Model transverse-momentum spectra of pions and kaons (solid lines) com- 
pared to the experimental PHENIX data (dots) [20] . 

the 3D transition temperature T$f is identified with the freeze-out temper- 
ature and the isotropization hypersurface is identified with the freeze-out 
hypersurface. 

Certainly, our method used to change from the purely transverse hy- 
drodynamic expansion to the isotropic system and further to freeze-out is 
very much simplified. One may think of several developments, for example, 
the transition criterion may be changed and/or the time evolution of the 
isotropic phase may be extended. Such developments are the subject of our 
current studies. 

5. The results of our numerical calculations with the initial distribu- 
tions corresponding to the non-central (c = 20-30%) Au+Au collisions at 
the beam energy ^snn = 200 GeV are shown in Figs. [2] - [H The values 
of the initial parameters used in the calculations are: Ti\ = 590 MeV, 
T^f = 170 MeV, and no = 0.43. We observe that our model reproduces 
well the transverse-momentum spectra of pions and kaons, Fig. [21 and also 
the pion+kaon elliptic flow, Fig. [31 both measured in Au+Au collisions by 
PHENIX [201 EU- In Fig- H we show tne model pionic HBT radii (solid 
lines) compared with the experimental STAR data (dots) [22]. The radii 
were calculated with the help of THERMINATOR which uses a two parti- 
cle method with Coulomb corrections [23J. The overall agreement with the 
data is quite good. The theoretical radii are by about 20% smaller than 
the experimental values, but their kx dependence is very well reproduced. 
Moreover, the model ratio i? out /i? s ide matches very well the experimental re- 
sult. Slightly smaller theoretical values indicate that the real isotropization 
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Fig. 3. Elliptic flow coefficient of pions and kaons (solid line) compared to the 
experimental PHENIX data (dashes) [21]. 



and freeze-out processes take longer time. This is, of course, an expected 
feature. 

For comparison, in Fig. 0] we also show the HBT radii obtained in 
the transverse-hydrodynamics model without the isotropization transition 
(dashed lines) [23]. We observe that the inclusion of the isotropization 
transition improves the kx dependence of R out and -R s id > which leads to 
the correct behavior of -R ou t/-R s ide- The main improvement, however, can 
be seen in the behavior of -Ri on g> which for small momenta is increased by 
a factor of 3. This means that the full 3D equilibration of the transversally 
thermalized system is necessary to reproduce the HBT radii. 

6. We conclude with the statement that the transverse- hydrodynamics 
model supplemented with the sudden isotropization transition and freeze- 
out is able to describe uniformly the soft pion and kaon data (transverse- 
momentum spectra, V2, and the HBT radii). This brings further evidence 
that the assumption of the very fast 3D equilibration of matter produced 
in the relativistic heavy-ion collisions may be relaxed. An earlier study 
showed that the RHIC data may be explained by the model where the 
standard hydrodynamic phase is preceded by the parton free streaming [25J . 
Our analysis of the HBT radii indicates, however, that the ultimate 3D 
equilibration is necessary to reproduce the experimentally measured value 

of -Rlong 

When we turn to heavier particles (protons, hyperons) we reproduce rea- 
sonably well their px spectra but find difficulties in the reproduction of their 
elliptic flow — the model results describing the elliptic flow of protons and 
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Fig. 4. Model results for the pionic HBT radii (solid lines) compared to the ex- 
perimental STAR data [22] (dots). The radii are plotted as the functions of the 
average transverse momentum of the pion pair kx- The model parameters are the 
same as in Figs. [2] and EH The dashed lines describe the results obtained without 
the isotropization transition |24) . 



hyperons are similar to the pion results, not exhibiting a distinct mass split- 
ting. Thus the model presented here does not explain the HBT puzzle under- 
stood as the consistent description of the soft observables for pions, kaons, 
and (at least) protons. On the other hand, the transverse-hydrodynamics 
model eliminates the assumption about early thermalization, makes a di- 
rect connection to color glass condensate, and describes consistently the 
pion+kaon sector with very economic tools. A further improvement of the 
model may be achieved by the matching with the energy-momentum tensor 
of the dissipative hydrodynamics. 

This investigation was partly supported by the Polish Ministry Of Sci- 
ence and Higher Education grant No. N202 034 32/0918. 
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